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ABSTRACT 

We measure and analyze the energy, momentum, and mass feedback efSciencies due to radia- 
tion from active galactic nuclei (AGN) in relatively large scale outflows (from ~ 0.01 to ~ 10 pc). 
Our measurements are based on the two-dimensional (axisymmetric) and time-dependent radiation- 
hydrodynamical simulations recently presented in Kurosawa & Proga. In that paper, we studied 
outflows from a slowly rotating (sub-Keplerian) infalling gas driven by the energy and pressure of the 
radiation emitted by the AGN. These simulations follow dynamics of gas under the influence of the 
gravity of the central lO'^ Mq black hole on scales from ^ 0.01 to ~ 10 pc. They self-consistently 
couple the accretion-luminosity with the mass inflow rate at the smallest radius (our proxy for the 
mass-accretion rate, Ma) . Over thirty simulations have been performed to investigate how the results 
depend on the gas density at the outer radius, po- A key feature of these simulations is that the 
radiation field and consequently the gas dynamics are axisymmetric, but not spherically symmetric. 
Therefore, the gas inflow and outflow can occur at the same time. We compare our Afa-po relation 
with that predicted by the Bondi accretion model. For high luminosities comparable to the Eddington 
limit, the power-law fit (A/a oc p^) to our models yields q w 0.5 instead of g = 1.0 which is predicted 
by the Bondi model. This difference is caused by the outflows which are important for the overall 
mass budget at high luminosities. The maximum momentum and mass feedback efficiencies found 
in our models are ^ lO^'^ and ~ 10^^, respectively. However, the outflows are much less important 
energetically: the thermal and kinetic powers in units of the radiative luminosity are ^ 10~^ and 
~ 10~^, respectively. In addition, the efficiencies do not increase monotonically with the accretion 
luminosity but rather peak around the Eddington limit beyond which a steady state disk-wind-like 
solution exists. Our energy feedback efficiencies are significantly lower than 0.05, which is required in 
some cosmological and galaxy merger simulations. The low feedback efficiencies found here could have 
significant implications on the mass growth of super massive black holes in the early universe. We 
stress however that we have not considered the innermost parts of the accretion and outflow where 
radiation and matter interact most strongly. The feedback from this region could have efficiencies 
significantly above the low values found here. 

Subject headings: accretion, accretion disks - galaxies: jets - galaxies: kinematics and dynamics — 
hydrodynamics - galaxies: evolution 



1. INTRODUCTION 

The central location of AGN in their host galaxies and 
the fact that they can produce a large amount of en- 
ergy imply that AGN can play a very important role in 
setting the physical conditions in their vicinity as well 
as on larger, galactic and even intergalactic scales (e.g., 
Igumenshchev et al. 1993; Ciotti & Ostriker 1997, 2001, 
2007; King 2003; Murray et al. 2005; Sazonov et al. 2005; 
Springel et al. 2005; Begelman & Nath 2005; Hopkins 
et al. 2005; Wang et al. 2006b; Thacker et al. 2006; 
Fabian et al. 2006, 2008; Pelupessy et al. 2007; Krohk 
2007; Merloni & Heinz 2008; Booth & Schaye 2009, and 
references therein) . There are many indications that sup- 
port this idea. For example, the presence of broad and 
narrow emission lines, broad and narrow absorption lines, 
in AGN spectra suggests that AGN continuum radiation 
affects the immediate environment of AGN (see Krolik 
1999 for an overview). In addition, the tight correlation 
between the mass (A/bh) of the central black hole (BH) 
in a galactic nucleus and the velocity dispersion a of the 
galaxy's bulge or spheroid, the so-called "A/bh — c" re- 
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lation (e.g., Ferrarese & Merritt 2000; Gebhardt et al. 
2000; Tremaine et al. 2002) can be explained by the feed- 
back between AGN and the infalling material from large 
distances. This feedback can quench both BH accretion 
and star formation in the galaxy when BH reaches a cer- 
tain mass. AGNs could provide such feedback because 
they are very powerful sources of energy and momentum 
(e.g., Silk & Rees 1998; Blandford 1999; Sazonov et al. 
2005; Fabian 1999; Fabian et al. 2002; King 2003; Scan- 
napieco & Oh 2004; Murray et al. 2005; Springel et al. 
2005; Di Matteo et al. 2005; Booth & Schaye 2009). 

AGN are powered by mass accretion onto a super mas- 
sive black hole (SMBH). To illustrate how the growth of 
SMBH can be self-regulated and how AGN feedback can 
be characterized, let us first express the radiation lumi- 
nosity due to accretion as 

La = erc'Ma, (1) 

where we invoke the simplest assumption such that the 
luminosity is proportional to the mass accretion rate 
(A/a) and a radiative (or the rest-mass conversion) effi- 
ciency (cr). Both er and A/a are uncertain. For example, 
ranges from ~ 10~^ in a standard, radiatively efficient 
thin disk to ~ 10~^^ for spherically symmetric accretion 
from a low density medium (e.g., Shakura & Sunyaev 
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1973; Shapiro 1973; Meszaros 1975; Soltan 1982; Yu & 
Tremaine 2002) while the mass accretion rate depends 
on poorly constrained physical conditions and geometry 
at large distances from the black hole. 

A common method to estimate is to adopt the ana- 
lytic formula by Bondi (1952) who considered spherically 
symmetric accretion from a non-rotating polytropic gas 
with uniform density p^o and sound speed Coo at infin- 
ity. Under these assumptions, a steady state solution to 
the equations of mass and momentum conservation exists 
with a mass accretion rate of 

A/b = A47rr|pooCoo, (2) 

where A is a dimensionless parameter that, for the New- 
tonian potential, depends only on the adiabatic index 
(cf. Bondi 1952; Shu 1992; Frank et al. 1992). The Bondi 
radius, re, is defined as 



where G is the gravitational constant and M is the mass 
of the accretor. 

One can quantify AGN feedback by measuring its ef- 
ficiency in affecting the flow of energy, momentum, and 
mass. In this work, we consider only the energy and mo- 
mentum carried out by matter. The total energy feed- 
back efficiency et is defined as the ratio between the ac- 
cretion luminosity of the system (Eq. [1] ) and the sum 
of the kinetic power (kinetic energy flux) and thermal 
energy power (thermal energy flux) Pth, i-e., 

et = (Pk + Ptb) /ia . (4) 

Similarly, the kinetic energy and thermal energy feed- 
back efficiencies (ek and eth) are defined as 

Ck = a/ia (5) 

and 

eth = Pth/ia , (6) 

respectively. From these definitions, it obviously follows 
that et = ek -I- eth- 

The momentum feedback efficiency (cp) is defined as 
the ratio of the total wind momentum p„ to the total 
radiation momentum {Lg_/c), i.e., 

ep = Pw/ (-^a/c) . (7) 

Lastly, the mass feedback efficiency e,,! is defined as the 
ratio of the mass-outflow rate at the outer boundary 
-^out to the mass-inflow rate at the inner boundary Min, 
i.e., 

em = Mout/Min . (8) 

The most advanced studies of feedback effects were 
carried out by Springel et al. (2005) (SDH05 hereafter), 
Di Matteo et al. (2005) and Booth & Schaye (2009) (BS09 
hereafter). These studies used computer simulations of 
merging galaxies in which they linked local and global 
processes. This was possible because they adopted rela- 
tively crude phenomenological realizations of star forma- 
tion, radiative cooling in a complex multi-phase medium, 
BH accretion and feedback, and because their spatial res- 
olution is larger than r^. In particular, they assumed val- 
ues of the above introduced efficiencies instead of directly 
computing them. A main result of these simulations is 



that the A/bh ~ cr relation can be reproduced remarkably 
well and that the relation is insensitive to the gas fraction 
in the model galaxies. In addition, the BH mass appears 
to be little affected by the details of star formation and 
supernova feedback. 

Begelman & Nath (2005) argued that the results from 
simulations of merging galaxies suggest that the feedback 
regulating BH accretion operates on local scales, compa- 
rable to tb or closer in, rather than solely on the global 
scales usually considered (see also Murray et al. 2005). 
They also presumed that the insensitivity to gas fraction 
occurs in the galaxy merger simulations because the gas 
mass is somehow "maximized" on the scales where the 
accretion rate is determined. Thus on these scales, BH 
feedback can be much more important than that due to 
stars. If this is correct then for state-of-the-art models, 
the key feedback processes represent "subgrid" physics. 
This is a limitation of current models of AGN feedback 
in large scale cosmological and galaxy merger simulations 
because they cannot be directly related to AGN physics. 

On the other hand, simulations that aim to provide in- 
sights to AGN physics do not include galaxy but rather 
focus on TB or even smaller scales (e.g., Proga 2007; 
Proga et al. 2008; Kurosawa & Proga 2008; Kurosawa 
& Proga 2009a, b). Thus they cannot be directly re- 
lated to AGN feedback on large scales. However, these 
smaller scale simulations can be used directly to measure 
the feedback efficiencies listed above. Consequently, they 
can be used to quantify the effects that are assumed or 
parametrized in large scale simulations. 

The goal of this paper is to present measurements of 
the mass accretion rate, and various feedback efficien- 
cies based on direct simulations of inflows and outflows 
in AGN, on sub-parsec and parsec scales performed by 
Kurosawa & Proga (2009a) (KP09 hereafter). In other 
words, we wish to determine if AGN can supply energy 
in the form and amount required by the cosmological and 
galaxy merger simulations. In this work, we do not at- 
tempt to provide a definitive answer to the problem of 
AGN feedback efficiency as our simulations do not in- 
clude the smallest scales (the black hole radius) and the 
very large scales (> 10 pc), and also do not include all 
physical processes operating in AGN (e.g., dust, mag- 
netic fields and star formation). Here, we simply report 
the AGN feedback efficiencies found in the simulations 
previously presented by KP09 who focused on relatively 
large scale inflow and outflow ( ^ 0.01 to ^ 10 pc). 

2. AGN MODELS IN CURRENT COSMOLOGICAL 
SIMULATIONS 

2.1. Mass Accretion Rates 

In recent cosmological and galaxy merger simulations 
(e.g., SDH05; Robertson et al. 2006; Sijacki et al. 2007; 
Khalatyan et al. 2008; Di Matteo et al. 2008; Johans- 
son et al. 2009; BS09), the actual physical process of the 
mass-accretion onto the BH is not explicitly modeled be- 
cause of a relatively poor resolution. Often, these simula- 
tions rely on a separate analytical model to describe the 
small scale physical processes. The unresolved accretion 
process is usually described by a Bondi-Hoyle-Little for- 
mulation (Hoyle & Lyttleton 1939; Bondi & Hoyle 1944; 
Bondi 1952). Here, we consider the case in which the ac- 
creting BH does not move with respect to the surround- 
ing gas. This reduces the process to a simpler Bondi 
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spherical accretion problem. Then, the mass-accretion 
rate can be written as Equation (2). The dimensionless 
constant A (in Eq. [2]) depends on the adiabatic index 7. 
For the gas with 7 = 5/3, A = 1/4 (see e.g., Frank et al. 
1992). The Bondi accretion formula relates the mass- 
accretion rate of a BH located at the center to the gas 
density and the sound speed (or equivalently the temper- 
ature) of the gas at a large scale. 

On the other hand, in the Bondi accretion prescription 
used by e.g., SDH05 and BS09, the mass-accretion rate 
is written as 



where p and Cg are the density and the sound speed esti- 
mated near the BH using the surrounding smoothed par- 
ticle hydrodynamics (SPH: Gingold & Monaghan 1977; 
Lucy 1977) gas particles. Note that the expression con- 
tains "the dimensionless parameter" a which is differ- 
ent from A in Equation (2). SDH05 introduced a pa- 
rameter to overcome the gap in the scale sizes between 
the numerical resolution and the Bondi accretion regime. 
In a typical cosmological or galaxy merger SPH simula- 
tion, the smoothing length (~ 10'^ pc) is much larger 
than the gravitational radius of influence or the Bondi 
radius, tb, which is 2 pc. If we assume the gas lo- 
cated at a large distance, heated by the AGN radiation, 
is "Comptonized" (T « 2 x 10^ K), the corresponding 
speed of sound (assuming 7 = 5/3) is relatively high 
(~ 500kms-i). 

SDH05, BS09 and others (e.g., Robertson et al. 2006; 
Sijacki et al. 2007; Khalatyan et al. 2008; Di Matteo 
et al. 2008; Johansson et al. 2009) find that a very large 
factor of a is required for low-mass BHs to grow their 
masses; hence, the problem is not strictly a Bondi ac- 
cretion problem. Most of the AGN feedback model in 
the cosmological simulations mentioned above assume a 
constant value of a = 100 (see also Table 2 in BS09) 
except for BS09 who allow a to depend on the value of 
local gas density. The assumption of a very large value 
of a becomes inadequate when the local gas density is 
higher than that required by the formation of the a cold 
interstellar gas phase, and when the cosmological simu- 
lation does resolve the Jean length and the Bondi radius 
(BS09). 

By noting these, BS09 abandoned the assumption of 
constant a in Equation (9), and introduced the following 
parametrization of a. 



1 for uh < njlj 

(nn/nn)'^ forn.n > 



(10) 



where nn and njj are the number density of hydrogen 
and the critical hydrogen number density above which 
the gas is expected to become multi-phase, and star for- 
mation is expected to begin via contraction of gas due to 
thermo-gravitational instability (cf. Schaye 2004; Schaye 
& Dalla Vccchia 2008). The critical density is chosen as 
nj^ = 0.1 cm i.e., the corresponding critical hydrogen 
density is pj^ = ^■'^ ^ 10~^^gcm~^. The best fit models 
of BS09 to some observations (e.g., the A/bh-c relation) 
gives (3 = 2.0. Note that the new parametrization of a 
in Equation (10) provides an additional density depen- 
dency of the mass-accretion rate in Equation (9). In the 



formulation of BS09, the mass-accretion rate steeply de- 
pends on the density of the surrounding gas i.e., M (x p'^ , 
for p > p"^ while the formulation of SDH05 and the orig- 
inal Bondi accretion model (Eq. [2]) always give a linear 
dependency i.e., M oc p. In most of the AGN accre- 
tion models in the cosmological simulations (e.g., SDH05; 
BS09), the highest mass-accretion rate is limited to the 
Eddington rate, i.e., 

47rGA'fBH'7ip 



(11) 



where 



Cr, (Tt and c are the proton mass, the ra- 



diative efficiency (the rest mass to radiation conversion 
efficiency), the Thomson cross-section and the speed of 
light, respectively. The Eddington ratio (F) is defined as 
the ratio of a system mass-accretion rate to the Edding- 
ton rate, i.e., F = Afa/A/Edd- 

Figure 1 illustrates how the mass-accretion rate de- 
pends on the density (p) in the models by Bondi (1952), 
SDH05 and BS09. The mass of the BH is assumed 
as A^BH = 10^ M0. In all three models, the speed 
of sound Cs is set to 520kms~^, which corresponds to 
that of Comptonized gas temperature T « 2 x 10^ K 
with the adiabatic index of gas 7 = 5/3. In the mod- 
ified Bondi accretion models of SDH05 and BS09, the 
mass-accretion rates are limited to the Eddington rate 
(Eq. [11]), and the radiative efficiency Cf is set to 0.1. 
The dimensionless parameter a ~ 100 is adopted for 
the model of SDH05, and f3 = 2.0 is adopted in the 
model of BS09. The figure shows that the mass-accretion 
rate of SDH05 is larger than the Bondi accretion rate 
for p ^ 10~^°gcm~'^. On the other hand, the mass- 
accretion rate by BS09 is similar to that of the Bondi 
accretion rate (off by a factor of A = 1/4 in Eq. [2]) in 
the low-density regime ( p ^ 10~^^ gcm~'^), but it is sig- 
nificantly larger than the Bondi accretion rate for the 
density range of 10"^^ ^ p ^ lO^^^gcm""^. The densi- 
ties above which the accretion proceeds at the Eddington 
rate are 9.5 x 10^^'^gcm^'^ and ~ 6.3 x lO^^^gcm^'^ 
for the models of SDH05 and BS09, respectively. How- 
ever, these values change depending on the adopted value 
of Cg. In the Bondi accretion model, the mass-accretion 
rates reaches the Eddington rate at much higher density 
(p ~ 10~^° g cm~'^). We note that adding a rotation to 
the gas can reduce the mass-accretion rate. As shown 
by Proga & Begelman (2003) and Ryu et al. (1995), the 
reduction of can be significant, i.e., by one order of 
magnitude or more, compared to Mb 

2.2. AGN Feedback Models 

In SDH05 and BS09 (also in many others, cf. Table 2 
in BS09), Er = 0.1 (e.g., Shakura & Sunyaev 1973; see 
also Soltan 1982) is assumed, and is fixed. A higher 
value of ej. (~ 0.2) can be achieved in an accretion 
model with a thin disk and a rapidly rotating BH (e.g., 
Thorne 1974). Recent observations suggest a wide range 
of e,.: 0.07 (Martfnez-Sansigre & Taylor 2009), 0.30-0.35 
(Wang et al. 2006a), 0.16 (Yu & Lu 2008), 0.15 (Elvis 
et al. 2002) and ^ 0.1 or ~ 0.2 (Yu & Tremaine 2002). 
On the other hand, Cao & Li (2008) find e^- is relatively 
low (-- 0.08) for A'/bh < 10® Mq and relatively high 
(<; 0.18) for A/bh <^ 10^ Mq. Since the exact mechanism 
of how the accretion luminosity of a BH couples to the 
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P (gcm"^) 



Fig. 1. — Comparison of the mass-accretion rates from the Bondi 
accretion model (Bondi 1952) (solid), the modified Bondi accretion 
model by SDH05 (dashed) and that by BS09 (dotted), as described 
in Eqs. (2), (9) and (10), respectively. The mass-accretion rates are 
computed as a function of density p at a radius ro which is much 
larger than the Bondi radius (re), i.e., ro 2> tb- In all models, 
the speed of sound Cs at ro is set to 520kms~^, which corresponds 
to that of Comptonized gas temperature T Ri 2 X 10^ K with the 
adiabatic index of gas 7 = 5/3. The mass of the BH is set to 
Mbh = 10* Mq. While a = 100 is adopted in the model of SDH05, 
/3 = 2.0 is adopted in the model of BS09. In the modified Bondi 
accretion models of SDH05 and BS09, the mass-accretion rates are 
limited to the Eddington rate (Eq. [11]). 

surrounding gas is not well known, SDH05 and BS09 sim- 
ply assume that La couples only thermally (and isotropi- 
cally) to the surrounding. Using Equation (1), the rate of 
energy deposition to the surrounding (the AGN feedback 
rate) in SDH05 is written as 

Ef ^EfL,^ etCrMBHc" (12) 

where et is the efficiency of the AGN energy deposition 
to the surrounding gas, and is a free parameter which is 
to be constrained by observations. BS09 find the models 
with et = 0.15 matches observations (e.g., the Magorrian 
relation and the Mbh — c relations) very well, and sim- 
ilarly SDH05 find et ~ 0.05 matches observations well 
(see also Di Matteo ct al. 2005). In the study of BS09, 
they find that the global BH number density at the cur- 
rent era (zero redshift) and the BH scaling relations are 
very sensitive to a choice of et, and the former is nearly 
inversely proportional to the value of et. 

3. OUR MODEL 

Our approach is to use physical two-dimensional (ax- 
isymmetric) and time-dependent hydrodynamical (HD) 
simulations of AGN flows to investigate the dependency 
of the BH mass-accretion rate on the surrounding gas 
density, and to find the AGN feedback efficiencies in con- 
verting the accretion luminosity into the outward fiuxcs 
of energy, momentum and mass. Here, we simply ana- 
lyze the simulations results previously published in KP09 
for this purpose. KPG9 used a modified version of the 
ZEUS-MP code (cf. Hayes et al. 2006) for their numer- 
ical simulations. In the following, we briefiy summarize 
their main model assumptions and results. 



In KP09, we studied axisymmetric hydrodynamical 
simulations of a slowly rotating gas that is under the 
influence of the gravity of a 10^ Mq black hole and is ir- 
radiated by a geometrically thin UV accretion disk and a 
spherical X-ray corona. We assumed that the AGN radi- 
ation is dominated by the disk radiation (95% of the total 
luminosity). Further, we account for the fact that the ra- 
diation from the disk depends on the polar angle 0, i.e., 
proportional to cos 9 due to the geometrical foreshorten- 
ing. We ran a set of simulations for various values of the 
gas density {po) at the outer radius of the computational 
domain, ro ~ 7 pc. After the initial transient stage, this 
density determines the key characteristics of our solu- 
tions such as the accretion luminosity and the outflow 
properties. We compute the accretion luminosity of a 
system based on the accretion-rate which is assumed to 
be equal to the mass-supply rate at the inner radius of 
the computational domain w 10~^ pc, (i.e., we used 
Eq. [1] where we assumed = 1/12 and Ma = Min[ri]). 
See § 5.3 for more detail on this assumption. 

For the models with high temperature gas at large radii 
and with high luminosities, we found a strong correlation 
between Mout and La (see Fig. 1 in KP09). The power 
law index describing the correlation is very similar to 
that for radiation-driven stellar and disk wind models 
(e.g.. Castor et al. 1975; Proga et al. 1998; Proga 1999). 
More surprisingly, for the models with high density at 
large radii, we found that steady state solutions with 
La exceeding the Eddington limit. The super-Eddington 
accretion proceeds in the equatorial region and is possible 
because the radiation flux from the disk is significantly 
reduced in the equatorial direction due to the geometrical 
foreshortening effect. 

In all simulations performed by KP09, an outfiow is 
driven from an infiow with sub-Keplerian rotation. For 
the models with high temperatures at large radii, the 
inflow occurs over a wide range of the polar angles, 
whereas the outflow occurs in a relatively narrow polar 
angle range (see the left panel in Fig. 2). However, for 
the super-Eddington cases with low temperature at large 
radii, the inflow persists only very close to the equatorial 
plane, resembling a thin accretion disk, while the outflow 
arises in a wide range of radii and polar angles (sec the 
right panel in Fig. 2). The geometry of this extreme 
inflow-outflow solution is very similar to a radiation- 
driven wind from a luminous Keplerian accretion disk 
(e.g.. Woods et al. 1996; Proga et al. 1998; Proga & 
Kallman 2002). For the cold super-Eddington solutions, 

Mout is only very weakly correlated. The weaker corre- 
lation is mainly caused by a mismatch between with the 
direction of escaping photons and the inflowing gas. In 
other words, the radiation is emitted mostly in the polar 
directions whereas the inflowing gas occurs mainly in the 
equatorial region. 

As it has been discussed and shown in the past, we 
find that self-consistently determined preheating/cooling 
from the quasar radiation can significantly reduce the 
rate at which the central BH is fed with matter. However, 
our results also emphasize a little- appreciated feature. 
Namely, quasar radiation does drive a non-spherical, 
multi-temperature and very dynamic fiow. 

In the following, we present the mass-accretion rates 
and various (energy, momentum and mass) AGN feed- 



5 



back efficiencies computed from tlic simulations. For tliis 
purpose, we use a subset of the models in KP09. Here, 
we concentrate on the models in which the outer bound- 
ary temperature is not fixed at a constant value, but it is 
self- consistently determined from the radiative and adi- 
abatic heating (Models 28-34 in KP09). Note that the 
flow solutions for these models are, in general, very sim- 
ilar to those with the fixed outer boundary temperature 
at 2 X 10^ K (the low temperature models, i.e.. Models 1-9 
in KP09). 

4. RESULTS 

We analyze the dependency of the mass-accretion rate 
on the gas density at a large distance from a BH and 
AGN feedback efficiencies based on the axisymmetric hy- 
drodynamical simulations presented in KP09. The main 
results of the models along with the input outer bound- 
ary density arc summarized in Tabic 1. 

4.1. Mass- Accretion Rates 
The mass-inflow rates at the inner boundary Mi^ (ri) 
and those at outer boundary Mi^ (ro) from the HD sim- 
ulations arc plotted as a function of the outer bound- 
ary density po in Figure 3 (see Tab. 1 for the numerical 
values). For a given value of po, Mn (ri) and Afin (ro) 
are not equal to each other, but rather Min (n) < 

Min (''o) because of an outflow. The lowest density model 
(Model 35) is an exception since no outflow is formed 
in this model. For the higher density models, an out- 
flow forms, and not all the material entering from the 
outer boundary reaches the inner boundary. A fraction 
of gas experiences a strong radiation pressure and radia- 
tive heating, and the direction of flow changes, forming 
an outflow. 

The figure also shows the the mass-accretion rates pre- 
dicted by the Bondi accretion model (Eq. [2]) and those 
computed from the formulations of SDH05 and BS09 
(Eqs. [9] and [10]). The outer radius is much smaller than 
that of a typical smoothing scale on a SPH cosmological 
simulation (~ 10'^ pc), and the outer density values used 
in our simulations are much larger than a typical local 
density at BH in the SPH simulations. In our simula- 
tions, the higher density at a 10 pc scale is required for 
a system to produce an outflow. For example, as we can 
see in Table 1, po must be greater than 2 x 10~^^ gcm~^, 
which corresponds to (tt-h ^ 1.2 x lO^cm"^), to form 
an outflow with our system setup. In the density range 
of the models considered here, the mass-accretion rates 
adopted by SDH05 and BS09 are limited by the Edding- 
ton rate (Eq. [11]); hence, the line is flat (cf. Fig. 1). Note 
that the radiative efficiency e,. = 1/12 instead of 0.1 is 
adopted for the Eddington rate (Eq. [11]) in the models of 
SDH05 and BS05 to be consistent with our simulations. 
This moves the Eddington rate only slightly upward. 

Our models include the effects of radiative heating and 
radiation force. Therefore, we do not, in general, expect 
our solution to reproduce an exactly same density de- 
pendency of the mass-inflow rate as that of the Bondi 
model. However, the figure shows the mass-inflow rates 
from our models are very similar to those of the Bondi 
rates, i.e., the rates are of the same order of magnitude. 
Interestingly our Mi^ (ri) and the Bondi mass- accretion 
rate matches around po = 4 x 10"^'' g cm~^ which coinci- 



dentally corresponds to P ~ 1. Since the accretion rates 
from SDH05 and BS09 are the Eddington rates (the rates 
corresponding to P = 1) in this density range, their lines 
also cross at the same point. 

The figure clearly shows that our models have a weaker 
dependency of the mass-inflow rates on the density than 
that of the Bondi accretion. The power-law fits of data 
points give the slope q = 0.52 (±0.01) for Min (n) and 
q = 0.56 (±0.02) for A/in (ro), which are indeed much 
smaller than that of the Bondi accretion model, i.e., 9 = 1 
(cf. Eq. [2]). Although not shown here, if we turn off 
the rotation, radiation force and radiative heating in our 
models, we obtain q « 1, as this is equivalent to the 
Bondi accretion problem (see also Proga & Begehnan 
2003; Janiuk et al. 2008). 

4.2. Feedback Efficiencies 

Next, we compute AGN feedback efficiencies in energy, 
momentum and mass using the simulation results, as de- 
fined in Equations (4)-(8). Since the models used here 
show some degree of variability (typically ~ 10% level, 
cf. Tab 1), the physical quantities used to compute the 
feedback efficiencies are based on the time averaged val- 
ues. 

4.2.1. Energy Feedback Efficiency 

Figure 4 shows the energy feedback efficiencies, et, Ck 
and eth computed based on our models (Table 1), as a 
function of the Eddington ratio (P). The numerical val- 
ues of the efficiencies are listed in Table 2. For systems 
with relatively low Eddington ratio (P^ 0.4), the ther- 
mal feedback efficiency is higher than the kinetic feed- 
back efficiency (eth > ek)- On the other hand, for sys- 
tems with relatively high Eddington ratio (P^ 0.6), the 
kinetic feedback dominates the thermal feedback by a 
factor of ~ 10 to - 100. The model with P = 0.2 does 
not form an outflow, indicating an approximate P value 
below which no outflow forms (with our system setup). 

The energy feedback efficiencies increase as P increases, 
but the efficiencies saturate for P 1. The total energy 
feedback efficiency peaks at P w 1 with et ~ 10~^. The 
flattening of the efficiencies for P ^ 1 is caused by the 
transition of the inflow-outfiow morphology to a "disk 
wind" like configuration (cf. Fig. 2) for the higher P mod- 
els (KP09). As briefiy mentioned in § 3, because of the 
mismatch between the direction in which most of the 
radiation escapes (in polar direction) and the direction 
in which the most of the accretion occurs in the system 
(the equatorial direction) , the radiatively driven outflows 
in the disk-wind-likc configuration cannot increase the 
outfiow efficiency by increasing the accretion luminosity 
or equivalently P. A similar behavior is found in the 
Mont (ro)- P relation of KP09 (sec their Fig. 7). 

4.2.2. Momentum Feedback Efficiency 

Figure 5 shows the momentum feedback efficiency ep 
as a function P. The numerical values of ep for each 
model are listed in Table 2. The dependency of Cp on P 
is similar to that of the energy feedback efficiency. For 
P ^ 1, ep increases as P increases, but it decreases as P 
increases for P ^ 1. No outfiow is formed for the mod- 
els with P = 0.2 and below. The cause of the peaking 
of ep at P « 1 (and the declining for P ^ 1) is again 
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Fig. 2. — Examples of the density and velocity maps from the set of simulations (without the temperature constrained at the outer 
boundary) performed by Kurosawa & Proga {2009a). The density (in logarithmic scale) is over-plotted with the directions of poloidal 
velocity as arrows for Models 36 {left panel), 29 (middle panel) and 34 {right panel). See Table 1 for the corresponding model numbers. 
The figures are placed in order of increasing density at the outer boundaries (po) from the left to right (cf. Table 1). Their corresponding 
accretion luminosities in the unit of the Eddington luminosity, i.e., the Eddington ratio, are F = 0.32,0.71 and 4.3, from the left to right 
panels. For the low density (po) and low accretion luminosity model (Model 36 in the left panel), the outflow is very narrow {0° <e < 30°), 
and the inflow is very wide (30° ^ 9 ^ 90°). As the density (po) and accretion luminosity increase (Model 29 in the middle panel), the 
outflow becomes wider (0° < 6 50°), whereas the inflow becomes narrower (50° ^ e ^ 90°). For the very high density and accretion 
model (Model 34 in the right panel), the outflow occurs over a very wide range of the polar angle (0° < ^ 85°), and the accretion region 
is now conflned to a thin equatorial wedge (the disk- wind- like solution). 



TABLE 1 
Model Summary 





Po 


T** 


Min in) 


F 


A^out {ro) 


Afin (ro) 


Model* 


(10-21 gcm-3) 


(10? K) 


(I025gs-1) 




(10^5 gs-l) 


(10^5 gs-i) 


35 


2 


0.02 


3.4(0.1)t 


0.20(0.01) 


0.0(0.0) 


3.4(0.1) 


36 


4 


0.07 


5.5(0.6) 


0.32(0.04) 


0.11 (0.02) 


5.4(0.2) 


28 


10 


0.14 


8.9(0.68) 


0.52(0.032) 


1.0(0.032) 


9.4(0.1) 


29 


20 


0.14 


12(0.58) 


0.71 (0.034) 


3.5(0.52) 


15(0.3) 


30 


40 


0.36 


18(1.1) 


1.1(0.1) 


7.1 (1.4) 


25(0.2) 


31 


80 


0.80 


25(2.6) 


1.4(0.2) 


10(2.3) 


35(2.1) 


32 


160 


0.98 


36 (4.9) 


2.1 (0.3) 


9.9(2.6) 


49(5.3) 


33 


320 


0.85 


52(2.6) 


3.1 (0.2) 


11 (0.48) 


63(1.1) 


34 


640 


1.30 


72(0.56) 


4.3(0.03) 


9.5(0.19) 


82(0.9) 



(*) The model numbers are identical to those in Kurosawa & Proga (2009a) except for Models 35 and 36 which are additional models 
presented here for the first time. 

(**) Self-consistently determined temperature at the outer boundary. 

(t) Values in brackets are the standard deviations of the time series values. 



TABLE 2 
Feedback Efficiencies 



Model 


F 


ek 








et 








em 


35 


0.20 

























36 


0.32 


1.4 X 10" 


8 


9.1 X 10" 


8 


1.0 X 10" 


7 


7.8 X 10" 


5 


0.02 


28 


0.52 


8.1 X 10- 


7 


4.6 X 10" 


-7 


1.3 X 10" 


6 


1.3 X 10" 


3 


0.11 


29 


0.71 


1.4 X 10" 


5 


1.0 X 10" 


-6 


1.5 X 10" 


5 


8.3 X 10" 


3 


0.29 


30 


1.1 


9.0 X 10" 


5 


1.8 X 10" 


-6 


9.2 X 10" 


5 


2.9 X 10" 


2 


0.39 


31 


1.4 


1.3 X 10" 


4 


9.0 X 10" 


-6 


1.4 X 10" 


4 


3.0 X 10" 


2 


0.40 


32 


2.1 


1.0 X 10" 


4 


2.0 X 10" 


-6 


1.0 X 10" 


4 


1.8 X 10" 


2 


0.28 


33 


3.1 


5.8 X 10" 


5 


1.7 X 10" 


-6 


6.0 X 10" 


5 


9.8 X 10" 


3 


0.21 


34 


4.3 


8.2 X 10" 


5 


1.5 X 10" 


-6 


8.4 X 10" 


5 


7.7 X 10" 


3 


0.13 
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Fig. 3. — Comparison of the mass-inflow rates found in the HD 
simulations (Tab. 1) with those predicted by the Bondi accretion 
model (Bondi 1952) {solid line) and with those adopted by SDH05 
and BS09 {dashed line). In the density range of the models consid- 
ered here, the mass-accretion rates adopted by SDH05 and BS09 
are hmited by the Eddington rate (Eq. [11]); hence, the hne is flat 
(cf. Fig. 1). The mass-inflow rates at the inner boundary {circles) 
and those at the outer boundary {squares) of the computational 
domain are shown as a function of the outer boundary density po- 
The mass-accretion from the HD simulations are very similar to 
the Bondi accretion rates, but the HD models has a less steeper 
dependency on the density. The Bondi mass-accretions rates and 
those of SDH05 and BS09 are computed for the gas with the Comp- 
tonized temperature T = 2 X 10^ K and with the adiabatic index 
7 = 5/3. 

due to the change in the inflow-outflow morphology to a 
disk- wind-like (KP09; see also Fig. 2). The momentum 
deposition of the photons becomes less efficient once the 
flow has the disk-wind-like configuration since a major 
fraction of the radiation escapes in the polar direction, 
but gas is not present in that direction since it has been 
already blown away by the strong radiation. The maxi- 
mum momentum feedback efficiency is tp w 0.03 which is 
about 2 orders of magnitude larger than the total energy 
feedback efficiency found in § 4.2.1. 

4.2.3. Mass Feedback Efficiency 

Figure 6 shows the mass feedback efficiency (em) plot- 
ted as a function F. The dependency of Cm on F is 
very similar to that of the momentum feedback efficiency 
(Fig. 5). The numerical values of em are listed in Table 2. 
For F < 1, the mass feedback efficiency em increases as 
F increases, but it starts to decreases slightly beyond 
F w 1. The efficiency peaks around F = 1 with the max- 
imum efficiency value ^ 0.4. In other words, about 40% 
of the total inflowing mass is redirected to an outflow. 
The turn-around of em values is caused by the transition 
of the outflow morphology to a disk-wind-like configura- 
tion (cf. Fig. 2) for the larger F models, as in the cases for 
the energy and momentum feedback efficiencies (§§ 4.2.1 
and 4.2.2). 

The accretion process in our model is fundamen- 
tally different from that of the Bondi accretion model 
and those adopted in the cosmological simulations (e.g., 
SDH05; BS09) because in our model an outflow and an 



Fig. 4. — The efficiencies of converting the BH accretion luminos- 
ity La to the rate of energy deposition to the surrounding gas are 
plotted as a function of the Eddington ratio (F). The panel shows 
the kinetic energy feedback efficiency ej^ {circles), the thermal en- 
ergy feedback efficiency {squares) and the total energy feedback 
efficiency et = ek + ^th {triangles), separately (see Eqs. [4], [5], and 
[6]). The maximum total energy feedback efficiency is ~ 10~^. For 
the models with relatively low Eddington ratio (F^ 0.4), the ther- 
mal feedback is more efficient than the kinetic feedback (etjj > cy^). 
For the models with relatively high Eddington ratio (F ^ 0.6), the 
kinetic feedback is more efficient than the thermal feedback by a 
factor of ^ 10 to ~ 100. The model with F = 0.2 does not form 
an outffow, and the vertical line {dashed) at F = 0.2 indicates an 
approximate F value below which no outffow forms. The flatten- 
ing of the efficiencies at F i?i 1 is caused by the transition of the 
inffow-outffow morphology to a "disk wind" like configuration for 
the larger F models (cf. Fig. 2). 
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Fig. 5. — The momentum feedback efficiency (ep), which is 
defined as the ratios of the total wind momentum (at the outer 
boundary) to the total radiation momentum (La/c), is plotted as 
a function of the Eddington ratio (F). The efficiency peaks at 
F Si 1 with ep ^ 0.03, and it decreases for larger F values. The 
decline of the curve beyond F fa 1 is caused by the change in the 
inffow-outffow morphology to a "disk wind" like configuration for 
the larger F models (cf. Fig. 2). The model with F = 0.2 does not 
form an outffow, and the vertical lino {dashed) at F = 0.2 indicates 
an approximate F value below which no outflow forms. 



w 

=§ 




0.01 - 



Fig. 6. — The mass feedback efficiency (em) plotted as a func- 
tion the Eddington ratio (F). The efficiency is defined as the 
ratio of the mass-outflow rate at the outer boundary A/out (to) to 
the mass-inflow rate at the inner boundary A/in C^i)i 

i.e., Cni — 

A'/out (^"0) /Afin (^i). The efficiency peaks around F = 1 with the 
maximum efficiency value ~ 0.4, i.e., 40% of the total inflowing 
mass is converted to the outflows. The turn-around of values 
is caused by the transition of the outflow morphology to a "disk 
wind" like configuration for the larger F models (cf. Fig. 2). The 
model with F = 0.2 docs not form an outflow, and the vertical line 
(dashed) at F = 0.2 indicates an approximate F value below which 
no outflow forms. 

inflow can simultaneously be formed while the Bondi ac- 
cretion model can only form an inflow. The model of 
SDH05 and others can form either an inflow or an accre- 
tion (but not both simultaneously). 

5. DISCUSSION 
5.1. Comparisons with Other Studies 

A main purpose of this study is to measure the ef- 
ficiency of the AGN feedback in our radiation hydro- 
dynamical simulations of accretion flows irradiated by 
AGN. Especially, we are interested in the energy feed- 
back efficiency (ef) as this is an important parameter 
to determine the BH growth rate and the BH number 
density at the current epoch (redshift zero) in the cos- 
mological simulations (e.g., SDH05; BS09). BS09 find 
that the BH growth rate is nearly inversely proportional 
to the value of ef , which is important for a self- regulation 
of a BH growth (e.g., Silk & Rees 1998; Fabian 1999). 

The approximate ratios of the maximum feedback effi- 
ciencies in energy, momentum and mass using the results 
from S 4.2 are: 



ep : et : ek : eth = 1000 : 100 : 1 : 1 : 0.1 



(13) 



where the total energy feedback efficiency et ^ 10~^. 
Compared to the mass feedback efficiency, the thermal 
energy feedback efficiency is 10* times smaller. Thus the 
coupling between the radiation and the thermal energy 
of the gas, on scales between 10^^ and a few parsecs, is 
relatively inefficient. 

The efficiency of the AGN energy deposition to the 
surrounding gas ef is usually treated as a free parame- 
ter and is somewhat constrained by observations, e.g., by 
fits to the Magorrian relation and the A/bh — c relations. 



SDH05, Robertson et al. (2006), Sijacki et al. (2007), 
Di Matteo et al. (2008) and Johansson et al. (2009) find 
that their numerical simulations with ef = 0.05 match 
observations. BS09, who adopted a slightly different BH 
accretion model (§ 2.1), find a larger value ef = 0.15 pro- 
duces a good match with observations. In their models, 
the energy is assumed to be released only in the form 
of thermal energy; hence, their ef is equivalent to our 
eth- Our thermal energy feedback efficiency eth is about 
5 X 10"^ times smaller than their ef . Even when we com- 
pare it with our total energy feedback efficiency et , their 
ef is still about 5 x 10^ times larger. 

The relatively small value of the kinetic energy feed- 
back efficiency ek found here (^ 10"''') is consistent with 
those found in Chclouche (2008) who analyzed the X- 
ray spectra {Chandra/ HETG) of five AGNs (type-I) with 
A/bh ~ 10^ M0. He found Ck ranges between 10~^ and 
10"'^. Our models also agree with the relatively low en- 
ergy and mass feedback efficiencies found by Krongold 
et al. (2007) and Stoh ct al. (2009) who studied of X-ray 
and UV absorbing outflows in Seyfert galaxies. On the 
other hand, the AGN evolution and the BH growth syn- 
thesis model of Merloni & Heinz (2008), combined with 
the observed star formation rate history of the universe, 
suggests that ek = (3 — 5) x 10"'^, which is about an order 
of magnitude larger than the value found in our models. 

While SDH05 and BS09 assumed that the energy feed- 
back efficiency is constant (independent of accretion lu- 
minosity), recent studies by Ciotti et al. (2009) and Shin 
et al. (2009), who performed one-dimcnsional hydrody- 
namical simulations of co-evolution of SMBH and ellipti- 
cal galaxy, include the dependency of the AGN feedback 
eSiciency (in a form of kinetic energy) on the accretion 
luminosity or more precisely on the Eddington ratio (F) 
of the system. Their parametrization of ek with F some- 
what resembles our result (Fig. 4), in a sense that ek 
monotonously increases for F ^ 1. Ciotti et al. (2009) 
find that models with varying ek arc in general more con- 
sistent with observations than those with a constant ek. 
Interestingly, ek values found by Shin et al. (2009) are 
only about 5 to 10 times larger^ than our values. This 
is a much better agreement than that with the models 
of SDH05 and BS09, as mentioned above. In summary, 
Ciotti et al. (2009) and Shin et al. (2009) concluded that 
the models need both radiation and mechanical feedback 
mechanisms included at the same time to account for ob- 
servations such as the ratio of the central BH mass to the 
stellar mass ratio (A/bh/A/*), and the X-ray luminosity 
of hot diffused gas. 

5.2. Causes of Low Feedback Efficiencies 

Relatively low feedback efficiencies found in our models 
(§ 4.2) are mainly caused by the combinations of the fol- 
lowing: (1) the axi-symmetric nature of the assumed un- 
derlying accretion luminosity (the disk radiation), (2) the 
axi-symmetric nature of the inflow-outflow geometry, and 
(3) the relatively low optical depth in the outflow. When 
the Eddington ratio (F) of the system is small, the inflow 
is nearly spherical. However, when the accretion rate is 
increased and F becomes close to 1, the flow becomes 

^ The feedback efficiency defined by Ciotti et al. (2009) and Shin 
et al. (2009) is equivalent to our X where fti 0.1 used in our 
model. 
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very aspherical (cf. Fig. 2). The inflow occurs only near 
the equator and the outflow occurs in the polar direc- 
tion. Despite a rather large value of T, the accretion can 
proceed in the equatorial region. This is because of the 
angular distribution of the radiation: in our model, the 
most of the radiation is emitted by a flat disk. Conse- 
quently, the disk radiation has cos 9 dependency due to 
the geometrical foreshortening, i.e., the flux peaks in the 
polar direction and decreases toward the equator. The 
Eddington ratio, F, is a global property of a system. The 
ratio {Q) between the radiation force (/r) and the grav- 
itational force (/g), i.e. Q = /r//g does not depend on 
the polar angle (6) for a spherically symmetric case; how- 
ever, it does depend on 9 for a case in which radiation is 
emitted from flatten source. In the spherical case, when 
F > 1, Q is always greater than 1 and an outflow forms 
in all directions {9). N.B. for the argument here, we con- 
sider the radiation force due to electron scattering only, 
but not including the force due to line processes (which 
is included in our simulations). For the disk radiation 
case, Q depends on 9, and Q can be less than 1 (inflow 
can occur) near the equator even if the global parame- 
ter F is greater than 1. In our simulations, Q > 1 in 
the polar direction (forming outflow) and Q < 1 in the 
equatorial direction (forming inflow) for all the models 
presented here, except for Model 35 which does not form 
an outflow. 

Once the inflow-outflow geometry becomes very as- 
pherical (i.e., equatorial) as in the middle and right pan- 
els of Figure 2, the coupling between the matter and 
radiation in the polar funnel is reduced. Most of the 
(high density) gas accretes near the equator where the 
radiation is weakest, whereas most of the radiation es- 
capes (without interacting with gas) in the polar direc- 
tions where there is very little accreting matter. The 
electron optical depth in the outflow in our models are 
rather small (tos < 1) even in the highest mass-accretion 
rate model. This mismatch between the preferred direc- 
tion of the disk radiation and the direction of accretion 
is a main cause of the relatively low feedback efficiencies 
found in our simulations. This occurs in a equatorial 
system, but not in a spherically symmetric system. 

In our models, the outflow is formed from initially in- 
flowing gas. In other words, some fraction of the in- 
flowing gas is turned around by the radiation force and 
becomes an outflow. The turn-around points or wind 
launching points can occur at a large distance from the 
center, i.e., near the outer boundary (cf. Fig. 2). When 
the wind is launched from a larger radius, the gas does 
not have time or space to accelerate to a terminal veloc- 
ity before escaping from the outer boundary. This may 
result in slightly smaller momentum and energy feedback 
efficiencies than those obtained from larger scale simula- 
tions. The eflaciencies reported here is strictly on ^ 0.01 
to ~ 10 pc scale. Finally, in the simulations presented 
by KP09, only gas dynamics and its microphysics are 
considered. Including the effect of dust would be very 
important in a larger scale simulation (> 10 pc). If dust 
dynamics is included in the model and the size of the 
computational domain is increased, the energy and mo- 
mentum feedback efficiencies could be larger than the 
values reported in this paper. 

5.3. Assumption on Mass Accretion Rate 



To compute the total luminosity of the system (Eq. 1), 
we need to have the information of the total mass ac- 
cretion rate (Ma) onto the BH. Since the inner radius 
(r; 0.01 pc) of the computational domain in the simula- 
tions in KP09 is much larger than the inner radius of the 
accretion disk (which is on the order of the Schwarzschild 
radius), we assumed that A/a is equal to the mass- inflow 
rate at the inner radius of the computational domain 
(Min [n])- In other words, we assumed that all the gas 
which crosses the inner boundary will eventually reaches 
the central black hole. We also assumed that no ma- 
terial can enter the computational domain from the in- 
ner boundary, whereas the energy and momentum in the 
form of the radiation can. These assumptions had to 
be made for numerical reasons. We simply do not know 
how much gas would enter from the inner boundary with- 
out explicitly modeling the accretion disk and its wind, 
consistently with the larger scale flow. Unfortunately, 
the current computer speeds and resources would not 
allow us to include such small scales and the large scale 
(~ 10 pc) flows at the same time in this type of numerical 
simulations. However, in reality, some gas would enter 
the computational domain from the inner boundary, per- 
haps due to a disk wind (e.g. Proga et al. 2000) or a jet 
produced in the immediate vicinity of the BH. We sim- 
ply do not know how much and in what form of material 
enter from the inner boundary; hence, we made the sim- 
plest assumption: no gas enters into the computational 
domain from the inner boundary. 

We recognize that the feedback from the innermost 
part of the accretion flow may be very important, and 
may influence the result of larger scale simulations such 
as the ones presented in this paper. Once again, we re- 
mind readers that the feedback efficiencies reported here 
are those of large scales (specifically ^ 0.01 to ~ 10 pc 
scales), but not in the AGN as a whole (including the 
accretion disk wind/jet). On the small scale, different 
forms of feedback mechanism (e.g. disk wind and jet) 
may be important, but this is beyond the scope of the 
current investigation. 

6. CONCLUSIONS 

We have presented and analyzed the AGN feedback ef- 
ficiencies of energy, momentum and mass based on our 
axisymmctric and time-dependent hydrodynamical sim- 
ulations (see Tab. 1) presented in Kurosawa & Proga 
(2009a). The simulations capture the radiation-driven 
outflows formed from a slowly rotating (sub-Keplerian) 
infalling gas under the influence of the gravity of the 
central SMBH. The accretion-luminosity and the outer 
boundary temperature arc sclf-consistcntly determined 
in these models. The radial range of the simulations 
spans from ^ 10~^ to ~ 10 pc. 

The dependency of the mass-accretion rate on the den- 
sity of the surrounding gas (or the outer boundary den- 
sity po) has been examined (Fig. 3). The result is com- 
pared with the Bondi mass-accretion rate (Eq. [2]), and 
with those adopted in the cosmological simulations of 
SDH05 and BS09 (see also Fig. 1). The density depen- 
dency of the mass-accretion rate in our models is some- 
what similar to that of the Bondi accretion model. For 
the density range of 10~^^ ^ Po lO^^^gcm"'^ (or cor- 
respondingly for 0.2 ^ F ^ 5), the differences between 
the two models are within a factor of 10. At F « 1 
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{po w 4 X 10^^° gcm~'^), the accretion rate of our model 
and that of the Bondi accretion model agree with each 
other. An important difference between the two mod- 
els is the steepness of the dependency on the density. 
The power-law fit (Ma Po) of o^'' models results yields 
q « 0.5 instead oi q = 1.0 which is predicted by the 
Bondi accretion model. This difference is due to out- 
flows in our model. We note that in this density range, 
the mass-accretion rates adopted by SDH05 and BS09 
have no density dependency because their accretion rates 
are limited by the Eddington rate. The accretion rates 
of SDH05 are artificially boosted up (a factor in Eq. [9]) 
by a factor of 100 in comparison with the Bondi mass- 
accretion rates. Consequently, their rates reach the Ed- 
dington limit at much smaller po than in our simulations 
(see Fig. 1). 

We find the energy feedback efficiency of our mod- 
els depends on the accretion luminosity of the system 
(Fig. 4). For F 1, the dependency is similar to the 
parametrization of the feedback efficiency adopted by 
Ciotti et al. (2009). Both kinetic and thermal energy 
efficiencies (ek and eth) peak at around F = 1. The max- 
imum efficiency values are Ck ~ 10~^ and eth ~ 10~^ re- 
spectively (see Table 2). For systems with relatively low 
Eddington ratio (F^ 0.4), the thermal feedback efficiency 
is higher than the kinetic feedback efficiency (eth > Ck)- 
On the other hand, for systems with relatively high Ed- 
dington ratio (F^ 0.6), the kinetic feedback dominates 
the thermal feedback by a factor of ~ 10 to ~ 100. 

The dependency of the momentum feedback efficiency 
ep on F is similar to that of the energy feedback effi- 
ciency (Fig. 5). The maximum efficiency is ep « 10^^ 
which is about ^ 100 times larger than that of the total 
energy feedback efficiency. The dependency of the mass 
feedback efficiency em on F is also similar to that of ep 
(Fig. 6). The maximum value of em found is ~ 0.4 at 
around F = 1, i.e., about 40% of the total mass that 
moves inward at large radii, does not reach the inner 
boundary of our computational domain, but rather is 
turned into outflows. 

Compared to the energy (thermal only) feedback effi- 
ciencies (ef — 0.05) required in the recent cosmological 
and galaxy mergers simulations (e.g., SDH05; Robertson 
et al. 2006, Sijacki et al. 2007, Di Matteo et al. 2008 
and Johansson et al. 2009), our thermal energy feed- 
back efficiency eth at the peak value is about 5 x 10^ 
times smaller than their ef . Our total and kinetic energy 
efficiencies are about 5 x 10^ times smaller than their 
values. These large discrepancies would suggest a few 
things. For example, our models are missing important 
elements. In particular, wc do not include effects of dust 
which could make the outflows stronger. In addition, 
we focus here on axisymmetric models which could differ 
from fully three-dimensional (3-D) models. Our prelimi- 
nary 3-D simulations show that in 3-D, the wind kinetic 
energy is smaller than in 2-D while the opposite is true 
for the thermal energy (Kurosawa & Proga 2009b). Al- 
though these changes are small (less than a factor of 2) in 
one of the cases studied by Kurosawa & Proga (2009b), 
they could be more significant in other cases (i.e., for the 
luminosity higher and lower than F = 0.6 assumed by 
Kurosawa & Proga 2009b). 



It is also possible that the AGN feedback may not be 
as effective as one might have had expected. Instead, 
other forms of feedback may be more significant than the 
AGN feedback via radiation on scales between 10~^ and 
a few parsecs. They include, the supernova feedback, the 
radiative feedback from star formation, the strong stel- 
lar wind from massive stars, and strong accretion disk 
winds or jets from AGN. The last two forms will intro- 
duce magnetic fields which may carry outward fluxes of 
energy and momentum. It is also possible that the AGN 
feedback efficiencies are indeed low and the AGNs take a 
long time to influence their environment. We note that 
in our models the AGNs do not shut off the mass supply 
completely even at very high luminosities. This indicates 
that the AGNs could operate on a very long time scale 
over which their impact on the environment can accumu- 
late, and eventually become significant. 

Finally, we conclude by noting that AGN feedback has 
two distinct modes: (1) radiation-driven (quasar) mode 
and (2) magnetic driven (radio jet) mode. In the cur- 
rent cosmological simulations, it is particularly difficult 
to deal with the latter, as it requires a full magnetohy- 
drodynamical (MHD) treatment, and we do not have an 
adequate resolution and more importantly the full un- 
derstanding of the jet mechanism. Nevertheless, Sijacki 
et al. (2007) extended the original implementation of the 
quasar mode feedback in SDH05 by injecting energy at 
random positions within a sphere centered around a BH 
to emulate the hot gas bubbles created by AGN jets. 
The radio mode could be very important in certain sit- 
uations, such as the cooling flows in clusters of galaxies 
where channeling the energy within a very narrow jet 
helps to transport the energy outside a galaxy. How- 
ever, this occurs in the low accretion rate regime, and it 
is subdominant in terms of the total BH mass growth. 
On the other hand, our simulations do not include the 
MHD treatment of the coUimatcd jet; hence, the focus 
of this paper is on the physical mechanism of the quasar 
mode feedback. Since the quasar mode is the dominant 
process for the total BH mass growth, our results of low 
feedback efficiencies would have significant implications 
on the mass growth of SMBHs in the early universe. In 
the near future, we plan to further investigate the impli- 
cations of our results by taking the boundary conditions 
of our simulations from full cosmological hydrodynamic 
simulations, and thereby making more direct connections 
with the physical conditions in a cosmological context. 
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